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ABSTRACT 

The X-ray cluster gas density distribution in hydrostatic equilibrium is computed 
from the universal density profile of the dark matter halo recently proposed by 
Navarro, Frenk and White (1996, 1997). If one assumes the isothermality, the resulting 
distribution is well approximated by the conventional /3-model. We predict the core 
radius r c , the /3-parameter, and the X-ray luminosity of clusters as a function of the 
temperature Tx of clusters in some representative cosmological models, and compare 
them with observations and results of numerical simulations. The predicted size of r c is 
a factor of (3 ~ 10) smaller than the average of observed values. If both the universal 
density profile and the hydrostatic equilibrium are reasonable approximation to the 
truth, then this suggests either that the previous X-ray observations systematically 
overestimate the core radius of gas densities in clusters of galaxies, or that some 
important physical mechanisms, which significantly increase the core radius, is still 
missing. 

Subject headings: cosmology: theory - dark matter - galaxies: clusters: general - 
X-rays: galaxies 
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1. Introduction 

While it has been recognized for a long time that the gas density profile of X-ray clusters is 
well approximated by the isothermal /3-model: 

n9{r) = [l + ( r /£)2]3/J/2' W 

the origin of this functional form is not yet accounted for. Navarro, Frenk, & White (1996, 1997; 
NFW96 and NFW97 hereafter) found that the virialized halos in numerical simulations are well 
fitted by the following universal profile: 



(r/r s )(l + r/r s 



pBM(r) = W1 C ° „/_ x 2 , (2) 



where p c o is the critical density of the universe at z = 0, and 

6 C (M) « 3 x 10 3 Q [l + Zf(M)f, (3) 



r(M) - !^M1 = J_( 3M V /3 (4 ) 
si ] ~ c(M) c(M)W^cpJ • [) 



In the above expressions, Oo is the density parameter at z = 0, A c (f2o>Ao) is the collapse factor 
in a spherical nonlinear model, 2f (M) is the average formation redshift of objects of mass M, and 
the concentration parameter c(M) is related to S C (M) as 



A c c 3 



3 ln(l + c) - c/(l + c) 



(5) 



NFW96 also mentioned that under the gravitational potential of the dark halo (0) , the gas in 
hydrostatic equilibrium reaches a profile quite similar to equation (|l]). Using a hydro-dynamical 
simulation, Eke, Navarro, & Frenk (1997, ENF hereafter) examined in more details the properties 
of X-ray clusters in a cold dark matter model with f^o = 0.3, Ao = 0.7, h = 0.7 and as = 1.05, 
where Ao is the dimensionless cosmological constant, h is the Hubble constant in units of 
100km-sec~ 1 -Mpc~ 1 , and as is the top-hat mass fluctuation amplitude at 8/i~ 1 Mpc. They also 
found that isothermal /3-model describes well the simulated cluster gas distribution. 

In the present Letter we derive an analytical expression for the gas density profile embedded 
in the universal dark matter halo (|2[) assuming hydrostatic equilibrium and isothermal distribution 
(§2). The resulting distribution is remarkably well fitted by equation (|l]). With this we are able to 
predict the core radius, the /3-parameter, and the X-ray luminosity of clusters as a function of the 
mass M or temperature Tx of clusters once the cosmological model is specified. Our predictions 
are compared with the results of ENF simulations and observations in §3. Finally §4 is devoted to 
discussion of further implications. 
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2. Hydrostatic equilibrium of gas and dark matter 



Consider an isothermal spherical gas cloud with temperature Tx, then its density distribution 
p g in hydrostatic equilibrium satisfies 



kT x dlnp g _ GM(r) 
fj,m p dr r 2 



(6) 



where \i and m p denote the mean molecular weight (we adopt 0.59 below) and the proton mass. 
If one neglects the gas and galaxy contributions to the gravitational mass in the right-hand side, 
then equation (0) yields that 



In 1 + 



M(r) = 4Tr5 c p c0 r 3 s 

and equation (Q) can be analytically integrated to give 

ln(l + r/r s 



r + r s 



Pgo exp 



with 



b(M) 



r/r s 

87rG(im p 5 c (M)p c or 



e -276/ 2 (1+r/rs) 27b/(2r/r fl)) 



27kT x 



(7) 



(8) 



0) 



The cluster gas temperature Tx is expected to be close to the virial temperature T v ; r (M) of 
the dark halo. In the profile (0), the latter is in fact dependent on the radius r: 



£T vir (r) = 7 



Gpm p M(r) 
3r 



(10) 



where 7 is a fudge factor of order unity which should be determined by the efficiency of the shock 
heating of the gas; ENF adopted 7 = 1.5 while Kitayama & Suto (1997) adopted 1.2 as their 
canonical value in the analysis of X-ray cluster number counts. If one substitutes equation (|io| ) 
into equation (||), one finds that 



6(r) 



2 r 
97 r s 



In 1 + 



r 



i -1 



r + r s 



(11) 



Throughout the present paper we assume that the cluster gas is isothermal; this is observationally 
true for most clusters, and theoretically expected due to the effect of the thermal conduction. In 
that case it is reasonable to adopt either T v ; r (r s ) or 2~vir(^vir) th.6 temperature Tx« This 
corresponds to b(r s ) w 1.15/7, or 6(r vir ) = 2c/(97)[ln(l + c) - c/(l + c)]" 1 w (1.3 ~ 2)/7 . 

The derived functional shape @ may appear completely different from equation (||) at a first 
glance. As is shown in Figure H, however, it is surprisingly well approximated by the isothermal 
/3-model profile: 



Pg{r) 



[l + (r/r CiCfr ) 2 ] 3 ^/ 2 ' 



(12) 
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with A(b) = -0.1786 + 0.982, r CjCfT = 0.22r s and f3 eS = 0.9b over 0.01r s < r < 10r s . Interestingly 
this implies that P e fi( r s) ~ 1/7 and /3 e ff( r vir)(l ~ 2)/7 if we take Tx = T v i r (r s ) and T v i r (r v i r ), 
respectively (also see Fig.Q b below). With 7 = 1 ~ 1.5, those values agree very well with typical 
observed ones, i.e., j3 = 1.2 ~ 0.6. This is why Waxman & Miralda-Escude (1995), NFW96 
and ENF were able to reproduce a gas distribution very similar to the /3-model in their specific 
models, and in fact the agreement should be quite generic and almost independent of the adopted 
parameters. 



3. Implications of the derived gas density profile 

Given the apparent success in describing the observed shape of the X-ray cluster gas profiles 
in the universal dark matter halo potential, it is reasonable to examine the predicted scales of the 
physical variables in further details. In fact, all the relevant quantities can be computed using 
Appendix A of NFW97 once the cosmological model is fully specified. 

Figure § plots the ratio of r v i r and r ce fj(= 0.22r s ) {panel a), the ratio of the virial 
temperatures at r s and at r v ; r {panel b), and r Cje ff {panel c), against the temperature Tx = T v ; r (r s ). 
In what follows we set 7 = 1.5 for definiteness, but it does not change the conclusions below. 
We adopt the shape parameter T to characterize the power spectrum, which is equivalent to 
£lohex.p[— f2&(l + V^h/flo)] in cold dark matter (CDM) models with the baryon density parameter 
Qb (Sugiyama 1995). The representative cosmological models which we consider include SCDM 
(too = 1, A = 0, r = 0.5, a s = 1.2), OCDM (fi = 0.3, A = 0, T = 0.25, a 8 = 1.0), and LCDM 
{n = 0.3, A = 0.7, r = 0.2, o- 8 = 1.0). 

Although r v j r /r c is around 50 on scales of rich clusters (Fig.||a), the corresponding 
temperatures T v ; r (r v i r ) and T vir {r s ) agree within a factor of two (Fig.|2|b). Therefore the assumption 
of isothermal clusters is acceptable, and one obtains reasonable values for /3 e g comparable to the 
observed one (§2) even if we adopt either T v ; r (r v i r ) or T v i r (r s ) as the gas temperature. In Figure 
||c, we also plot the r c and Tx relations of observed clusters from the Einstein (Jones & Forman 
1984; David et al. 1993), and EXOSAT (Edge & Stewart 1991a,b; see also Kitayama & Suto 
1996) observations. Clearly the predicted r Cje g is a factor of (3 ~ 10) smaller than the average of 
observed values. 

Equation (|6|) does not involve the amplitude of gas density, the central density p g o should be 
determined by some other consideration. For this purpose, we use 

-^gas ('Vir ) /gas ^ b 



M(r vir ) O 

where / gas is the gas mass fraction of the total baryon in the cluster. Observationally / gas is 
typically 0.7 ~ 0.9. With equations (0) and (§), equation (13) reduces to 



(13) 



Zgas^fePcO^c o 27f>/2 V ' 



fin 



-e- 



ln(l + c) 



1 + c 



f x 2 {l + x) 27b ' 2x dx 
Jo 



(14) 
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This can be translated to the central number density of electron n e o assuming the primordial 
abundance of hydrogen and helium (X = 0.76), which is plotted in Figure |3|a. For definiteness, 
we set b = 0.7/0.9 so as to reproduce a typical shape of observed X-ray clusters, and adopt 
£lb = 0.015/i~ 2 and / gas = 0.8. As is clear from Figure ||a, the predicted n e o is larger than typical 
observed values. This is a direct consequence of the small r c in this model given the total gas mass 
ratio in the cluster at r v ; r (Fig.||c). 

Once n e o is specified, it is straightforward to compute the X-ray bolometric luminosity: 

f'vir 

Lx,bo\{Tx) = 47r / r 2 dr a{Tx)n 2 e (r) 
Jo 

= 4ira(T x )n 2 e0 r 3 s e- 27b j° x 2 {\ + x) 27b/x dx, (15) 

Jo 

where a(Tx) denotes the bolometric emissivity for which we include line emissions in addition 
to the bremsstrahlung (Masai 1984, Kitayama, Sasaki & Suto 1998). The results are plotted in 
Figure |||b, together with the observed £x,boi ~ Tx relation (David et al. 1993; Ebeling et al. 1996; 
Ponman et al. 1996). Since -Lx,boi predicted from equation (|l5|) scales as h~ 3 unlike the observed 
ones (oc h~ 2 ), we divide our predictions by h assumed in each model in Figure |||b. Dash-dotted 
line indicates the best-fit to the observed relation 

£ X ,boi = 2.9 x 10 44 /i~ 2 ( —J erg sec" 1 , (16) 

(Kitayama & Suto 1997). The shallower slope of the predicted Lx,boi ~ Tx relation (~ 2) than the 
observed one (~ 3.4) reflects a self-similar nature of the evolution in this model (Kaiser 1986). 
In fact our model predictions reproduce the results of Kitayama &: Suto (1997; their Fig.l) and 
of ENF (their Fig. 15), apart from the small difference due to the different choice of adopted 
parameters. 

4. Discussion and conclusions 

One of the most important consequences of the universal density profile is that all the 
physical quantities characterizing the gas density profile (r c , r v i r , /3 e fr, ^x.boh ^eo) are computed as 
functions of the total halo mass M, or almost equivalently of the corresponding gas temperature 
Tx, once underlying cosmological parameters (Qq, Ao, T, h, and erg) are fully specified. This can 
be regarded as a significant theoretical improvement of the understanding in the physical origin of 
the conventional /3-model. 

On the other hand, if both the universal density profile and the hydrostatic equilibrium are 
reasonable approximation to the truth, then our result indicates either that the previous X-ray 
observations systematically overestimate the core radius of clusters of galaxies, or that we neglect 
some unknown important physical mechanisms which significantly increases the core radius. 
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Incidentally r v ; r in equation (||) would be further divided by 1 + Zf if Z{ really corresponds to the 
formation redshift of the entire cluster, and this would make r c even smaller. 

Figure |3] shows the mass profile of the universal dark matter halo. Also plotted are 
the estimates on the basis of the j3- model fitting (eqs.|jl2| and §). The mass inferred from 
the gas density in this way underestimates the true value for r < r c . In turn, if the X-ray 
observations overestimate the core radius of clusters of galaxies, they would underestimate the 
total gravitational mass of the clusters. 

ENF obtained that r vir = 2.1/i~ 1 Mpc, r s = O^/i^Mpc, r c = O.l/i^Mpc, and 
n e Q = 2 x 10 -2 /i~ 2 cm~ 3 by averaging over the ten most massive clusters at z = 0. Their 
r v i r , r s and r c are systematically larger than, but yet are consistent (within a factor of two) with, 
our LCDM predictions (Figs. ^ and ^); this might be partly ascribed to the insufficient number of 
particles to resolve the small core, and partly to the fact that they include clusters not completely 
in equilibrium. Fukushige & Makino (1997) and more recently Moore et al. (1997) argued that 
their simulations with much higher spatial resolution show the steeper inner density profile like 
oc r~ L4 which should predict the core radius much smaller than what we discussed above. 

We thank Yipeng Jing for providing numerical routines to compute 8 C (M, z) in various 
cosmologies. We are grateful to Tetsu Kitayama and Noriko Yamasaki for useful discussions, and 
to Tatsushi Suginohara for comments on the manuscript. This research was supported in part 
by the Grants-in-Aid for the Center-of-Excellence (COE) Research of the Ministry of Education, 
Science, Sports and Culture of Japan (07CE2002) to RESCEU (Research Center for the Early 
Universe), University of Tokyo, Japan 
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Fig. 1. — Gas density profile (solid lines) expected from the universal density profile of dark matter 
halo (dashed line) for b = 1/3, 1, and 5/3. For comparison, the best-fit /3-models with f3 = 0.96 
and r c = 0.22r s are plotted in dotted lines. 
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Fig. 2. — Predicted properties of gas density distribution as functions of the X-ray cluster gas 
temperature Tx = T v ; r (r s ). (a) Ratio of r v ; ir and r c (« 0.22r s ); (b) ratio of the virial temperatures 
at r s and at r v ; r ; (c) predicted sizes of the effective core radius compared with the observed cluster 
data. Solid, dotted and dashed lines indicate SCDM, OCDM, and LCDM models, respectively. 
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Fig. 3. — Central electron number density and bolometric X-ray luminosity plotted against the gas 
temperature for Q b = 0.015/i" 2 , / gas = 0.8 and b = 0.78 (/3 cfr = 0.7). 
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Fig. 4. — Mass profile of the universal dark matter halo (solid line) compared with the estimates 
on the basis of the /3-model fitting. Dotted line indicates the estimate for the best-fit /3-model 
(/3 = 0.96 and r c = 0.22r s ), while dashed line corresponds to = 0.96 and r c = r s - 



